TRACING KAM TORI IN PRESYMPLECTIC DYNAMICAL 

SYSTEMS 

HASSAN NAJAFI ALISHAH AND RAFAEL DE LA LLAVE 

Abstract. We present a KAM theorem for presymplectic dynamical 

systems. The theorem has a " a posteriori " format. We show that given 

fvj a Diophantine frequency u} and a family of presymplectic mappings, if we 

y—^ find an embedded torus which is approximately invariant with rotation 

f^ u) such that the torus and the family of mappings satisfy some explicit 

fSJ non-degeneracy condition, then we can find an embedded torus and a 

h^,^ value of the parameter close to to the original ones so that the torus 

jrt is invariant under the map associated to the value of the parameter. 

^H Furthermore, we show that the dimension of the parameter space is 

^H reduced if we assume that the systems are exact. 



1. Introduction 



O 

cn 

Cd Presymplectic structures (constant rank, closed 2-forms) arise naturally 

S i„ the ..udy of .e.ene..te La^a„.ia„ and Ha.„ilto„ia„ .elanical systems 

with constrains, in time dependent Hamiltonian systems and in control the- 
> ory. (see, e.g., (6}[8l[l2l[l3l[20l[2ll[27] ) . 

l/S Given a presymplectic form Jl G 17^ (M), a vector field X G X(M) is said 

^) to be a Hamiltonian vector field associated with a function H G C'^{M) if: 

cn 

O ix^ = dH. 

(N 

IJ Due to the degeneracy of 17, there can be different functions H associated 

k>< with X, not differing by a constant. The corresponding flow </)^ : M — )• M is 

^ a 1-parameter group of presymplectic diffeomorphisms: ((/)^)*ri = il. Hence, 

the dynamics of such systems leave the presymplectic structure invariant. 

For other situations where presymplectic dynamics occur see, e.g., (2]. 
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Our aim is to state and prove a KAM type theorem for presymplectic 
dynamical systems, which extends the results of iTj for the case of symplectic 
diff eomor phisms . 

Our main result can be stated as follows. Let T" = ]R"/Z" be the n- 
dimensional torus. We consider the presymplectic manifold 

(1.1) M := T*T'^ X T"" , 

with an exact presymplectic form 0, of rank 2d, whose kernel coincides 
with the T"'-direction. See Remark ll.2l and Remark l3.5l for more observations 
about the kernel. 

One says that K : T^" — )• M is an invariant torus of a diffeomorphism 
f : M ^ M with frequency oj € M"+"' if: 

f{K{e))-K{e + cj) = 0, V0eT"+'^. 

When the left hand side is non-zero, but small enough (in some smooth 
norm that will be made explicit later), one says that / has an approximate 
invariant torus K with frequency uj (this will be made precise later). Our 
main theorem can be stated in rough terms as follows: 



Theorem 1.1. Let f\: M ^ M, where M is as in (1.1), be an analytic, 
non- degenerate in the sense to be defined later, {2d + n) -parametric family 
of presymplectic diffeomorphisms such that /o has an approximate invariant 
torus Kq, satisfying a non- degeneracy condition, with frequency uj satisfy- 
ing a Diophantine condition. Then there exists a diffeomorphism fx^ in 
this family, where Aqo is close to 0, which has an invariant torus K^o with 
frequency oj and which is "close" to the initial torus Kq. 

The precise version of the theorem will be stated below in Section [3j We 
need to formulate precisely the non-degeneracy conditions and to make more 
precise the definition of "close" which requires introducing norms. Note that 
we will not require the system to be neither nearly integrable, nor to be writ- 
ten in action-angle variables. Indeed, just like in |7J for the symplectic case, 
the fact that the dynamics of the system preserve the presymplectic struc- 
ture implies that the KAM tori are automatically approximately reducible. 
This leads to an approximate solution of the linearized equations without 
transformation theory. Moreover, the reducing transformation is given ex- 
plicitly in term of the approximately translated torus, which form the basis 
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of an efficient numerical algorithm (an explicit description of this algorithm 
can be found in [16|). We will discuss all this in more detail in the Section 
[3] below. 

The proof of our main theorem follows an approach similar to the one 
developed in JTJ for the symplectic case. The presymplectic case however 
has a few peculiarities due to the degeneracy of the 2-form. For some of the 
most routine calculations, we will just refer to this paper. The quasi-Newton 



method used here (and in |7|) is of the type introduced by Moser in 23,24 



We note that the approach is not based on transformation theory, which 
seems problematic in the case of presymplectic mappings since generating 



functions are not as straightforward as in the symplectic case (see |4lp, 10 



for studies of the theory of canonical transformations) and the Lie transform 
method is hampered by the fact that there are several Hamiltonians that give 
the same vector field. The approach is based on deriving a parameterization 
equation and applying corrections additively. The presymplectic geometry 
leads to cancellations that reduce a Newton step to the constant coefficients 
cohomology equations customary in KAM theory. We also note that the 
same cancellations lead to very effective numerical algorithms. 

We will also prove a fiux-type vanishing lemma for exact presymplectic 
diffeomorphisms. Roughly speaking, we will show that the average of the 
translation is zero in the directions other than the ones tangent to torus in 
the basis. Note that in the directions tangent to torus the averaging does 
not need to vanish. This also shows the need for considering a parametric 
family of diffeomorphisms, rather than just a single diffeomorphism. 

There are two possible ways to extend the results, we have stated here 
for the maps, to flows. One is to use the local uniqueness we state in the 



theorem 3.7 and the other one is to proceed with automatic reducibility for 
flows as in the symplectic case which is done in |7]. Both of them require 
some technical details and it is work in progress to do so. 

As a final note, it should be remarked that the results of this paper are 
not applicable to the dynamics associated with general Poisson structures. 
For regular Poisson structures, which have an underlying regular symplectic 
foliation, there are cohomological obstructions to find a compatible presym- 
plectic structure, see [30]. Even when these obstructions vanish (e.g., lo- 
cally around invariant tori), so that one can find a compatible presymplectic 
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structure, Poisson diffeomorphisms do not coincide with presymplectic dif- 
feomorphisms, and these two kinds of diffeomorphisms have quite distinct 
properties. The KAM theory for Poisson manifolds has been developed in 



26 . In section 8l we will compare the Poisson and presymplectic cases. 



Remark 1.2. An important well known fact about presymplectic forms is 
that the kernel of Q is an integrable distribution. 
We recall that the kernel of a form is 

Ker(O) = {Z\iz{n) = 0} = {Z\n{Z,X) = VX} 

Note that, for a general 2-form fi and any three vector fields X, Y, Z, we 
have: 

dn{x, Y, z) = x{n{Y, z)) - Y{n{x, z)) + z{9.{x, y)) 

-^{[X,Y]Z) + n{[X,ZiY])-^{[Y,ZlX) 
If dVi = and X, Y are in the kernel of $7, for any Z, we have ^^([X, Y],Z) = 
0. 

Hence, ii X,Y e Ker(J7), [X,Y] € Ker(O). This shows that the distribu- 
tion given by the kernel can be integrated to a manifold. 

Of course, in a torus, it could well happen that the leaves integrating the 
kernel are not compact (e.g. they could be an irrational foliation). 

1.0.1. Some exam^ples and com,m,ents. One example to keep in mind could 
be the three dimensional torus endowed with a presymplectic form Vt = 
d^i A d^2- Clearly, the kernel is given by the level sets of 'I'l, ^2- 

A more complicated example on T^ is $7 = d^i A 7 where 7 is a closed 
but not exact form. In this case, the kernel can be an irrational foliation. 

Another example related to the previous ones is the study of quasi- 
periodically perturbed Hamiltonian systems H{x,ujt). These can be made 
autonomous by adding an extra variable E T"^ that satisfies -^6 = lj. The 
phase space is now supplemented by a factor T . The symplectic form in 
the phase space becomes a presymplectic form in the extended phase space 
having T'^ in the kernel. Even this elementary example was considered as 
covered by the KAM theory of symplectic systems at the time of writing 

EZl- 

The theory of presymplectic manifolds was developed (e.g. in [15]) to give 



a geometric framework to the Dirac theory of constrained systems, 118, 19 



There are many physically interesting examples of constrained systems to 
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which the present theory appHes. Notably, besides the examples in 18 , 19 



the papers 12, 13 contain a very concrete example of a relativistic system 
of spinning particles which is close to integrable. 

The paper 16] shows how the Pontriaguin maximum principle for optimal 
trajectories can be formulated using presymplectic systems. If we consider 
a mechanical system with KAM tori and subject it to a control indexed by 
enough parameters, the results in this paper give a condition which ensures 
that the one adjust parameters to maintain the quasi-periodic motion. It 
would be interesting to study in detail concrete models, specially because 
the methods of this paper are well suited for numerical implementations. 

Note that, in contrast with symplectic manifolds, presymplectic manifolds 
may be odd dimensional. Hence, it is clear that an extension of the symplec- 
tic theory to presymplectic systems will require significant modifications. A 
general theory of perturbations of quasi-periodic motions independent of 
geometric structures was undertaken in [25] . Note, however that the per- 



turbations by presymplectic systems do not satisfy the assumptions of 25 
because the normal eigenvalues of a torus do not change when we change 
the parameters. Of course, in the presymplectic case the origin of this de- 
generacy is the preservation of the geometric structure, which also helps by 
eliminating some of the perturbing terms. This is the interpretation of the 
geometric identities used in this paper: they kill several dangerous pertur- 
bation terms. 

2. Preliminaries and Motivation 

In this section we will fix some notations and state a few preliminary 
results. Along the way, we will also justify the assumptions that will appear 
later in our main result. 

As stated in the introduction, we consider M = T"' x M'^ x T" equipped 
with a constant rank exact presymplectic structure, i.e., an exact 2-form 
Q G Q,'^{M), such that its kernel is: 

N ■.= Kein = {{u,{0,0,z)) £TM \ ueM,z€R"'}. 

The exactness assumption places restrictions on the presymplectic form 
but for applications to Hamiltonian dynamical systems with constrains or 
degenerate Lagrangian systems this is not too restrictive. For these systems 
the phase space is often obtained by restriction to a submanifold where the 
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2-form is the pullback of the canonical symplectic structure on the cotangent 
bundle (see [ISi]) or some other exact symplectic form (see |20]). 

Let V = {{u, [x, y, 0)) G TM\u G M, (x, y) G M^d} ^^ ^i^^^ j,j^ = v ®N, 
and denote by vr : TM — )• V the canonical projection on V . For each 
z G TM, we have the linear isomorphism J{u) : T^M — t- T^M defined by: 

(2.1) nuiC,r]) = {^,J{u)7]), C,7]eTuM 
where 

and (•, •) denotes the standard Euclidean inner product on M^^". The 
skew-symmetry of il implies that J^ = —J. 

We will be using the following norms. If x = (xi, ..., Xd+n) £ W^~^^ we set: 

|x| := max \xj\. 

j=l,..,d+n 

For an analytic function g on a complex domain B we denote by |g|c™,B its 
C™-norm: 

|5|c™,B := sup sup\D''g{z)\, 
0<\k\j,<m zeB 

where |% := |/i| + ... + \ld+n\- 

We will be looking for real analytic invariant tori which extend holomor- 

phically to a small strip in the complex space. More precisely, let Up denote 

the complex strip of width p > 0: 

Up = {9€ C'^+"/Z'^+" : |Im(p)| < p}, 

and introduce the following family of maps. 

Definition 2.1. The space ("Pp, ||-||p) consists of functions K : Up ^ M 
which are one periodic in all their arguments, real analytic on the interior 
of Up and continuous on the closure of Up. We endow this space with the 
norm 

(2.2) ||i^||p:= sup|if(^)|, 

eaUp 

which makes it into a Banach space. 

We will also use the same notations for functions taking values in vector 
spaces or in matrices. 

Some well known results about the spaces above are the Cauchy bounds 
below (a consequence of Cauchy's integral representation of the derivative). 
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For < 6 < p, we have: 

(2.3) \\D^K\\p^5<Cj6-^\\K\\p 

Like in all other KAM type results we will have to deal with small divisors. 
For that we set: 

Definition 2.2. Given 7 > and a > d + n, we will denote by D{'y, a) the 
set of frequency vectors cj G R'^"*"" satisfying the Diophantine condition: 

(2.4) \l ■ LU - m\ > -f\l\^'' V/ GZ'^+"\{0},mGZ 

The aim of this paper is to find invariant tori of a given frequency u for a 
m-parametric family of presymplectic diffeomorphism fx, defined as follows: 

Definition 2.3. A m-parametric family of presymplectic difFeomor- 
phisms fx is a function 

f : M X B ^ M, 5 C M'", 

such that for each x £ M the map f{x,-) is of class C^ and for each A G -B 
the map fx '■= /(•, A) is a real analytic presymplectic diffeomorphism. 

We will introduce an algorithm to solve the equation 

(2.5) fx{K{e))-K{e + u;) = 0, uj G D{-f,a), 

given that one knows an approximate solution Kq{6) for the diffeomorphism 
/aq, where, with out loss of generality, we will set Aq = 0. In other word, we 
know that 

(2.6) fx,(.Ko{9))-Ko{e + uj) = eo{e), 



where the error term eo(^) has small enough norm. Equation (2.5) will be 
solved by Newton method where at each step we have infinitesimal equations 
given by 

dfx{K,i9)) 



(2.7) Dfx^iKiie))A,{e) - Aiie + uj) + 



£i = -ei{9). 



dX 
The approximate invariant tori and the geometry of the problem will lead 



us to a change of variables that will reduce (2.7) to a simpler equation with 
constant coefficients that can be solved by the following result of Riissmann 
which will also be useful in some other proofs. 
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Proposition 2.4 ((l6[[28]). Let uj G D{a,^) and assume that h : T"'+" -)> 
■^2d+n j^g analytic on Up and has zero average, avg(/i) = 0. Then for all 
< 6 < p , the difference equation 

(2.8) v{e) - v{9 + uj) = h{e) 

has a unique zero average solution v : T^" — )• M^ "*"" which is analytic in 
Up^g. Moreover, this solution satisfies the following estimate: 

(2.9) \\v\\p-5<co^-^6-''\\h\\p, 
where cq is a constant depending on n and a. 

2.1. Lagrangian properties of invariant tori. A first, very important, 
consequence of the Diopliantine condition on uj and the exactness of the 
presymplectic form is that invariant tori are actually Lagrangian tori: 



Lemma 2.5. If K{9) & Vp is a solution of (2.5) then K*0, is identically 
zero. 



Proof. Since K{9) satisfies (2.5) and / is presymplectic we have 

K*n = (KoT^yn, 

where T^{9) = 9 + oo. Moreover, since uj is rationally independent, rotations 
on the torus are ergodic and this implies that K*Q. is constant. If we write 



K*Q. in matrix form, exactly as we did for Vt in (2.1) we have 

(2.10) i^*f](C,r?) = (e,L(%) e,^GTe(T^+") 

where L[9) is actually constant. It remains to show that L{9) = 0. 
The 2-form Q. is exact, so we can write O = da where 

a{u) = a{u)du, a{u) = (ai(u), ..., a2d+n{u)y ■ 

Then we find that 

d+n 

{K*a) = ^Cj{9)d9^ 

where the components Cj have the following expression 

Cy9) = DK{9)a{K{9))j. 
This implies L{9) = DC{9y - DC{9). But now: 

(2.11) s.yg{DC{9)) := / DC{9)d9 = 0, 

which shows that: 

avg(L(0)) = 0. 
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But L{9) being constant, we conclude that L{9) = 0, i.e., K*Q = 0. D 



Following simple lemma extends the result of the Lemma 2.5 to approxi- 
mate invariant tori: 

Lemma 2.6. Let fo '■ M ^- M be a presymplectic analytic diffeomor- 
phism and let K G Vp be an approximate invariant torus with frequency 
uj e D(7,(t); 

(2.12) foiKi6))-Kie + oj) = eie). 

and assume that Jq extends holomorphically to some complex neighborhood 
of the image of Up under K: 

Br = {ze €2^^+" : sup \z - K{9)\ < r}. 
e&Up 



Then there exist a constant C > 0, depending on n, o", p, ||Z)i^||p, |/o|(7i g^ 
and \J\c^fi^, such that for Q < 5 < ^ 

(2.13) ||i||p-2a<C7-ij-("+i)||e||p 



where L is the matrix representing the puHback form K*il. (see 2.10). 



Proof. Let g := L — L o T^^. Then, we note that g is the expression in 
coordinates of 

K*n - T*K*n = K*f^n - T*K*n 

Hence, when K is exactly invariant g = 0. One can also easily show that 
llsll < ll-^e|| See [7| for more details. 



Using Proposition 2.4 one obtains that: 

ll-^llp-25 < Co7"-^(5"'^||5||p_5. 

One can bound the norm of g in exactly the same way as in the symplectic 
case, which can be found in [7J, to obtain the result. D 

2.2. Automatic reducibility near invariant tori. In this subsection we 



will assume that K{9) is an invariant torus of /, i.e., a solution of (2.5) 



When one starts instead with an approximate invariant torus Kq{9) of /, 



i.e., a solution of (2.6), the results of this subsection do not hold anymore. 
However, we will see in the next sections that we have versions of these 
results which hold in the approximate case and which will allow us to perform 
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the Newton method and conclude the existence of an invariant torus. For 
K{9) £ Vp let us decompose its Jacobian in the form 

(2.14) DK{9) = iX{e),Zi9)) 

where X(9),Z{9) are the first d and last n columns of DK(9). Also, for 
every vector in TM = V (B N, we will use the subscripts V and N for the 



first and second projections in each factor. Assume that K{9) solves (2.5) 
and that there exists a d x d-matrix valued function N{9) such that 

(2.15) Ni9)iXli9)-Xvi9))=ld, 



where X{9) is as in (2.14) This non-degeneracy assumption will turn out to 
be one of the ingredients to solve (2.5) approximately. Also, selll 

(2.16) Yv{e) := Xv{9)N{9) and Y{9) := ^^^^^ 

Then the following matrix will provide us the change of variable needed to 

reduce the linearized equations ( |2.7| ) to a simple form: 

Xvie) J-\K{9))Y{9) Zv{9) 
Xn{9) Zn{9) 



(2.17) 



M[9) 



where X,Z and Y are defined in (2.14) and (2.16) respectively. The non- 



degeneracy assumption (2.15), together with the fact that K{9) is Lagrangian 



(Lemma 2.5), show that: 

(2.18) ^K{e){X{9), J-\K{9))Y{9)) = h 

(2.19) nK(^e){X{9),X{9)) = Q 

(2.20) ^Kie){X{9),Z{9)) = Q 

Therefore, X{9), J^^{K{9))Y{9) and Z{9) do not form a presymplectic basis 
along the torus K{9), because neither nK(e){J'^{K{9))Y{9), J-\K{9))Y{9)) 
nor Q,j^(^0-^{J^^{K{9))Y{9), Z{9)) have to be zero, but they do provide a ba- 
sis where ri takes a rather simple form. Moreover, as the following lemma 



shows, they transform the linearized equations (2.7) into a simpler form: 
Lemma 2.7. The set {X{9), J-'^{K{9))Y{9), Z{9)} is a basis provided the 



matrix 
(2.21) 

V{9) = 



Id 

-/, -YH9)J-\K{9))Y{9) {J-\K{9))Y{9)yj{K{9))Zv{9) 
Xn{0) Zn{9) 



We will often abuse notation and will use Y{6) to denote both Yv{0) and Y{6). 
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is invertihle. In this case, we have: 
(2.22) 

Df{K{e)).{X{e), Z{e)) ={X{e + a;), Z{e + a;)), 
Df{K{9))J-\K{e))Y{9) =X{e + uj)Si{e) + J-\K{9 + u}))Y{e + a;) Id + 
(2.23) +Z{e + uj)A{9), 

where A[9) and Si{9) are matrices satisfying: M 
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(2.24) 

Proof. Let 
(2.25) 



Df{K{9)) ■ M{9)) = M{e + w) • 



h 


Si{9) 1 





Id 





A{e) In \ 



Q{9) :-- 



xl{e)j{K{e)) 

{J-\K{9))Y{9)YJ{K{e)) 







J-r>. 



The expression ( |2.17[ ) for M and relations ( |2.18[ ), ( |2.19[ ) and ( |2.20[ ), give: 
(2.26) 



Q(^)-M(0) 



/rf 

-/, -yT(0)j-i(i;c(0))y(0) {J-\K{e))Y{9))U{K{9))Zv{9) 

XN{e) z^(e) 



which shows that {Xv{9),J-^{K{9))Y{e)} is a basis for V := tt{TM) and 
Q{6) is invertible. Using this fact one can write 

Zv{e) = afie)XUe) + h^{9)J-\K{9))Y\ (/ = 1, . . . , d, A; = 1, . . . ,n). 

Pairing both sides with Xy{9) via the presymplectic form Q, it foUows from 



( |2.18D , ( |2.19D and ( |2.20D that: 

(2.27) 6f„(^) = f^(X(°(0),Z^(^))-af(0)(7(X|«(^),xU0)) 



0. 



In general, we have no control on Q.{J^^[K{9))Y{6),Zv{9)), it means we 
have no control on the af{9), but the assumption that V{9) := Q{6) ■ M{9) 
is non-degenerate guarantees that {X{9), J^^{K{9))Y{6), Z{9)} is a basis. 
Assume from now on that V{9), and hence M{9), is invertible. Since / 



is presymplectic and f{K{9)) = K{9 + u}), is follows from (2.18), (2.19) and 



We emphasize that identity (2.241 holds only when we have an invariant torus. In 



Corollary 4.2 we will prove that for approximately invariant tori (2.241 holds up to an 



error which can be bounded by the error in the invariance equation. 
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( |2.20D that: 

Df{Kie)).ixie), zie)) ={x{e + u), z{e + a;)), 

Df{K{9))J-\K{e))Y{e) =X{e + u)Si{9) + J-\K{9 + uj))Y{e + 6^) Id + 

+ Z{9 + uj)A{e), 

for some matrices Si{9) and ^(^). This shows that relations ( |2.23 ) hold. 
Moving the term J^^{K{9 + u}))Y{9 + u) Id to the left side of the second 
equation we obtain that: 

(2.28) 
A{e) = me + u) [Df{K{9))J-\K{9))Y{9) - J-\K{9 + uj))Y{9 + oj)] , 

where T3(0) is the last row in the matrix: 

Ti{9) 



(2.29) 



M-\9) 



T2{9) 



Finally, moving the term Z{9+uj)A{9) to the left hand side and pairing both 
sides with J~^{K{9 + iLi))Y{9 + w), via the presymplectic form il, together 



with (2.18), gives: 



(2.30) Si{9) = [Yv{9 + uj)y 0] [Df{K{9))J-\K{9))Y{9)- 

J-\K{e + uj))Y{9 + uj)- Z{e + uj)A{9)] . 

D 

Remark 2.8. A straightforward calculation shows that V^^{9) takes the 
following the form: 



(2.31) 



We will need this fact later. 



V{{ V{, V,-, 
Id_ 0_ 0^ 

^31 ^32 ^33 



3. Main results 

In this section we will give precise statements of our results. The discus- 
sion in the previous section motivates introducing the following: 

Definition 3.1. We will say that K{9) G Pp is a non-degenerate torus 
if 

(i) There exists a d x d-matrix valued function N{9) such that 

(3.1) N{9){Xy{9)y.{Xv[9))=^^ 
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(ii) the matrix V{6), which defined in (2.21), is invertible. 



where Xv{0) and V{9) are defined in (2.14) and (2.21) respectively. 



Remark 3.2. In the symplectic case, the matrix V{9) is always non-degene- 
rate when fi is exact. When il. is not exact, then one also needs to assume 
that V{9) is invertible in order to perform the Newton iteration successfully. 
In the presymplectic case, even when fi is exact, we need to assume that 
V{9) is invertible. Also, we do not know how to proceed with the algorithm 
presented here if one gives up on exactness of 0,. However, one may still be 
able to proceed with this algorithm in some special problems where the form 
is non-exact. Dealing with KAM theory for non exact symplectic forms is 
a deep problem largely unexplored, see [29] for remarks on the problem of 
non-exact forms. 

Definition 3.3. A pair {fx,K{9)) is non-degenerate at A = Aq if fx is a 

(2d -|- re)— parameter family of presymplectic diffeomorphisms, K{9) € Vp is 
a non-degenerate torus, and the average of the {2d -|- re) x {2d + re) matrix 

(3.2) A{9):=V-H9)Q{9)(^ {K{9))] 

A=Ao / 



has rank 2d + re, where V'(^) is defined by (2.21 ). 




Remark 3.4. As indicated in the outline above, the role of the non-degeneracy 
assumption is that we can transform the Newton equation to a constant co- 
efficient equation (up to a small error). 

Note that the condition is an open condition, so that, if the initial er- 
ror is small enough, the iterative process does not leave the region where 
Definition 13.31 holds. 



Remark 3.5. Note that the Definition 3.3 involves the presymplectic form. 



In geometric terms, the assumption in the non-degeneracy condition we 
need is that the tangent space to the torus, its symplectically conjugated 
space and the space of tangents along the parameter in the form span the 
whole space. That is, that the direction of moving along the parameter 
space, compensates the kernel of the presymplectic form. We have formu- 
lated this assumption using the fact that the kernel corresponds to some 
coordinates of the space due to the integrability of the kernel which comes 
for free from the closeness assumption of the presymplectic form. 
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We can now state the main theorem of the this paper: 

Theorem 3.6. Let uj G D('y,a), let fx be a 2d + n-parametric family of 
analytic presymplectic diffeomorphisms and let Kq £ 'Ppo- Assume that: 

(HI) The pair [fx^Ko) is non- degenerate at X = Xq. 

(H2) The family f\ can he holomorphically extended to some complex 
neighborhood of the image of Up under K: 

Br = {zeC: sup \z - K{e)\ < r] 

such that \f\\c'^^Br < ^^^- 

If 

eo{9):=fxo{Ko{9))-Ko{e + u;), 

then there exists constant c > 0, depending on a, n, d, po, r, I/aqIc^.b^' 
\\DKq\\p^, ||A''o||po; 11^ Ia=Ao (-f^o)llpo o.'^d |avg(Ao)~^| such that if < 
(5o < max(l, ^) and 

(3.3) lleollpo < min{7^5o'^,rc7^5o''l|eo||po} 

then there exists a mapping Koo G Ppo_6o-o and a vector Aoo £ M^d+n g^^j^. 
fying 

(3.4) fx^ oK^ = K^oT^ 
Moreover, the following inequalities hold: 

(3-5) ll-f^oo - -fi^ollpo-6.50 < -7 ^Q '^lleollpo 

(3.6) lAool < ^7''5o"'"l|eo||po 

Sketch of the proof. More details of the proof wih be given later, in Sec- 
tions 4][5][6, but it will be useful to start with a brief overview that can serve 
as a road map. 

We will use a modified Newton method of the type introduced by Moser 



in [23,124,31 . The procedure goes as follows. Starting with 

(3.7) G{Ko,Q) := fo{Ko{e)) - Ko{e + uj) = eo{e), 

we look for an approximate solution for the corresponding linearized equa- 
tion 



(3.8) DG{K,M(i,,,^e),e,) 
dfx{Ko{e)) 



dX 



£0 + DfoiKoie))Aoie) - Ao(0 + uj) = -eo(0). 

A=0 



KAM FOR PRESYMPLECTIC MAPS 15 

The left hand side of this equation By an approximate solution we mean up 
to a quadratic error, i.e., a solution Ao{6) such that: 

\\DG{Ko,0)\iA„ie),eo) + eo\\po-5o < co7"''Jo"^'"^'^l|eo||^o 
where 5o,co are constants to be determined later. 

Having the solution (Ao(0),eo) a better approximating torus for the map 
/ai, where Ai = Ao + eo, is defined as 

Ki{9) = Ko{e) + Ao{9) 

and it will be shown that {Ki{9), /a^) is a non-degenerate pair. Furthermore, 
setting 

we find that 

In other words, for the new torus the error has decreased quadratically. 
Iterating this procedure, we will see that the sequence 

(i^o,Ao),(i^i,Ai),...,(K„,A„),... 

of approximate solutions of ( [2.5^ , obtained by applying the iterative pro- 
cedure, converges to a solution (Kac^oo)- One has to be careful with the 
domain Up which decreases in each iteration (the reason is because we can 
bound the correction applied at one step only in a domain slightly smaller 
than the domain of the original function). This loss of domain can be ar- 
ranged in a way that, in the limit, one does not end up with an empty do- 
main. This choice of decreasing domains so that there is some domain that 



remains is very standard in KAM theory since the first papers 114 , 23 , 24 



See 16,31 for a pedagogical exposition. D 



3.1. Local uniqueness. Notice that if K^o is a solution of (3.4) then for 
every (/? G T x T" the map Koo{6 + ^) is also a solution. For this reason, we 
will consider K(9) and K{9) := K{9 + ip) to be equivalent. By uniqueness 
of solutions, we will mean uniqueness up to this equivalence relation. The 



following result gives uniqueness of solutions of (3.4): 

Theorem 3.7. Let oo G D{j,a) and assume that Ki and K2 are two non- 
degenerate tori in Vp solving 

(3.9) fx{K{9)) -K {9 + u;) = 0, 
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such that Ki{Up) C Br and K2{Up) C Br- Furthermore, assume that the 
matrix 

U.-avgl ^^^^ 



where Si{9),A{9) are defined by (2.30) and (2.28) has rank d. Then there 
exists a constant c > depending on a, n, d, p, r, I/aIc^.b^' ll-^-^illp? Il-^illp 
and \@\ such that if 

(3.10) \\Ki-K2\\p<cj^6^'', 

where 6 = §, then there exists an initial phase r G T x T" such that in Up/2 
one has: 

KioTr = K2 

The proof of this result is given in section [7] 

3.2. A vanishing lemma. We end this section with one geometric result. 
Recall that a diffeomorphism / : M — >• M is called exact presymplectic 
if at the level of de Rham cohomology one has: 

(3.11) [/*a-a]=0 

where a is a primitive of the presymplectic form: = da. When M is not 
compact, one must use compactly supported de Rham cohomology. Clearly, 
the time-1 map of a Hamiltonian vector field is exact. Moreover, using the 
flux homomorphism (see [I]), one can show that an exact presymplectic 
diffeomorphism which is close enough to the identity is the time-1 map of a 
(time-dependent) Hamiltonian vector field. 

We now generalize to exact presymplectic diffeomorphisms the Vanishing 
Lemma of J9], valid for exact symplectic diffeomorphisms, and which allows 
one to have some control on the size of the parameter A. Due to the the 
presence of kernel, our Vanishing Lemma has a slightly different nature (and 
statement) than [Oj Lemma 4.9]. 

We will assume that we are in the situation described in the statement of 



Theorem 3.6, where /o is exact. In order to simplify the notation we write 
K{6) instead of Koo{0) and A instead of Aqo- Let fx := fx — /o and define 
the averagq^ 

(3.12) fl:= I fx{K{d)) de £ R2'^+" 

Jjd+n 



In the sequel, we will not distinguish between a map with values in f''+" x R'' and a 
lift with values in R2''+". 
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If we express the vector /Z in the basis {X{6), J^^{K{9))Y{9), Z{6)}, we 
obtain the ^-dependent components (/ii(0), ..., fi2d+n{d)), in other word 

fi=[fi,{e),...,fi2d+nm[x{e) j-\K{e))Y{e) z{e)]. 

We have 

Lemma 3.8 (Vanishing Lemma). If fo : M ^ M is an exact presymplectic 
diffeomorphism, then 

(3.13) / fik{0) de = 0, {k = d+l,...,2d). 

Jfd+n 

Proof. We fix the following notations 

Gi = (^1, •••,^j-l,^i+li •••i^d+n) S T 

U)i = (wi, ...,(^j_i,a;j+i, ...,t0dj^n) £ I^ 
and we let Uj^e- : T — )• T'^^" be the path given by; 

'^iA^'n) = (Oi, ...,9i^i,r],ei+i, ...,9d+n)- 
Also, we consider the two-cell B^ ^ : [0, 1] x S^ — ^ ]g2d+n (^gflned by: 

(3.14) B^^^^{C,rj) ■.= Koa^^^^_^^X^) - (/2) o a,^^^^.^(r/)e 
We will compute the integral 



in two distinct ways: 

(1) The boundary of -B^g. is the difference between the two paths K o 
a- g , ^. and {K oT^^ — fl) o a-^, so by Stokes's theorem we conclude 

(3.15) n= a- a 

Since {f\{K{9))—fl) has average zero and satisfies all hypothesizes of Propo- 
there exists an analytic function v : T "*"" — t- M^ "''" such that 

v{e)-v{e + uj) = fxoK-fL. 



sition 



2.4 



This, together with the exactness of /o implies that: 

/ a = a = a 



[ /o«+/ 

J Kou. % J voa . 



a 



a + a. 

Koa. n J voa. s —voa. 



z,w,+^?' 
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Hence, we see that: 






a. 



By a simple change of variable, we see that if we integrate over the torus 
qpd+n-i ^]^g right-hand side of the previous equation vanishes, so we can 
conclude that 



(3.16) 



dOi 



Q = Q. 



(2) Next we compute the integral of O over B. ^ explicitly as follows: 



n 



B. 



1 /-l 



JO 



^B^^Xi,v){diB0^{i,r]),dr,B0^{i,vi))did7] 



Since ^ is a constant vector, by (3.14) and (2.14) we have for i = 1, . . . , d: 



dvB^ 



K o a ■ a , - = Xj o a ■ a I - 



9^Big^, = -ifl)oa.j 



ifii+uji' 



So from the partial presymplectic basis relations (2.18), (2.19) and (2.20) 
we conclude that: 



(3.17) 



VL= fid+i o cr ... (r/) dry, (i = 1 . . . , d). 



Now, (3.16) and (3.17) together show that: 



d9i fid+ioa.Q^^Xri) dri = {i = l...,d), 
and this yields the result. 



n 



Remarks 3.9. The following remarks illustrate the relevance of the Vanishing 
Lemma: 

• The Vanishing Lemma concerns invariant tori. It can be extended to 
the approximate case, as it is done in [9] for the symplectic case, and 
assuming that the whole family fx is exact presymplectic, it leads to 
a bound on the parameter, which shows that in every step the value 
of the parameter decreases with the error term. This can be useful 
in numerical schemes for finding invariant tori. 

• In dimension 2, a volume preserving diffeomorphism of S^ x R is the 
same as (pre) symplectic diffeomorphism. In this case, as shown by 



the proof above, the integral (3.13) is the oriented area between a 



circle and its image by the map, as shown in Figure [T] This clearly 
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shows that the vanishing of (3.13) is an obstruction for the existence 
of invariant tori (see also [16|). 




Figure 1. Vanishing Lemma 

• Recall that we can think of our presymplectic manifold M as T*T'^ x 
T". In our Vanishing Lemma we only control the averages in the 
directions normal to T . It is easy to give simple examples of maps 
satisfying all the assumptions and such that the averages in other 
directions are non-zero. 

4. Estimates for the linearized equation 



The sketch of the proof of Theorem 3.6, given in the previous section 



relied on finding an approximate solution of the linearized equation (3.8) 



assuming that one has an approximate solution Kq of (3.7). In this section, 
we explained how this can be done. 

The first claim is that the set {X{e), J-^{K{e))Y{e), Z{e)] is stih a basis 
for Txg(e)M if the error term is small enough. Note that now, due to the 



error term, equation (2.26) becomes 



(4.1) 
where 



Qie)-M{0) = v{e) + R{9), 



R{e) :-- 



xi{e)j{e)Xv{e) o xi{e)j{e)Zv{e) 





If we now use that Kq{6) is approximately Lagrangian, i.e., if we apply 



Lemma 2.6 we see that we can control the reminder R(9): 
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Lemma 4.1. Assume the hypotheses of Lemma \2. 6\ hold. Then there exits 
a constant C3 depending on d, n, p, |/a|ci.B; I-^Ic^B; II-^IIp; ^^^^ \\DKo\\p 
such that for every < 5 < § we have 



^-'•i?||p-25<C37-'5-('^+^)| 



eol 



p- 



We conclude that: 



Corollary 4.2. Assume the hypotheses of Lemma 2.6 hold. Ifeo{6) satisfies 

1 



C37-'r('^+i)||eo||p<-, 



(4.2) 

then M is invertible and 

M-\d) = v-\e)Q[e) + MM, 

where 

(4.3) M,{e) = -[hd+n + v-\e)R{e)]-^v-\e)R{e)v{e)R{e). 

Moreover 

(4.4) ||Me||p_25<C47"''^"^"+'^||eo||p, 

where C4 is a constant which depends on the same parameters as C3. 

Proof. A simple application of the Neumman series. See 1^. 

We are ready to apply our change of variables. Before that we remark 
that, since /aq is presymplectic, we have 

FiiO) 



D 



(4.5) 



DfxoiK{0)) 





F2i9) F^iO) 



where Fi{9) is a symplectic linear map from V = n{Tp^tQ\M) into itself. 

Lemma 4.3. Let Kq{6) € Vp solves 

fx,{Ko{e))-Ko{e + uj) = eo{e) 
and that {f\,K{6)) is non- degenerate at X = Xq in the sense of definition 



3.3 Ifeo{e) satisfies (4.2), then the change of variable Ao{0) = M{e)i{e) 



transforms equation (3.8) to 
(4.6) 



h s{e) 
h 

A{e) In 



+ B{e)\m-ao + o^) 



dfx 



v-He)Qi9)eo{e) - A(0)eo - M,{e)eo{e) - M^m"^ 



)£o, 



A=Ao 
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where 

' S2{e) 




E{e) := {Dieo{e),Ei{e),D2eo{9)) 



B{e) ■.= M-\e + oj)E{e) 



E,{9) := Dh,{Ko{e))J-\Ko{e))Y{6) - X{e + io)Si{e)+ 

- J~HKo{e))Y{e + u;)- Z{6 + Lo)A{e) 

S2{e) ■= VI3 • {F2{e)j-\Ko{e))Y{e) - Xn{9 + oj)Si{e) - z^ie + uj)a{9)) 
s{e):=Si{e) + S2{e), 



andA{9), Me{6) and Si{9) are definedby (3.2), (4.3) and (2.30) respectively. 
Moreover, we have the estimates: 



(4.7) 

Me 

(4.8) 



|MeeoL-25<C47"'r('^+i)| 



eollp 



d{fx o Ko) 



dX 



£0 



A=0 



151 



p-25 
p~25 






d{fx o Ko 



dX 



A=0 



|eo|||eo||p 



eol 



where C4 is the same as in (4.4) and C5 is another constant which depends 
on the same parameters . 



Proof. The form of the transformed equations follows from substituting the 
change of variable and elementary computations. 



To prove the estimates (4.7) and (4.8), we note that (4.7) follows imme- 



diately from (4.4), so it only remains to prove (^4.8|. First note that for 



the first term in the definition of B{9) i.e. M~^{6 + co)E{9), the Cauchy 
provide bounds for L'ieo(^) and D2eo{9) in terms of the error. This enables 
us to bound M~^{9 + uj) {Dieo{9), D2eo{9)) by the error term. Calculating 
bounds for M^^{9 + uj)Ei{9) is more subtle. By the definition of A[9) and 
the fact that 



it follows that 



T3{9 + u)X{9 + u}) = 



n{9 + u;)Ei{9) = 0. 



Therefore: 



M-^{9 + uj)Ei{9) 



Ti{9 + oj) 
T2{9 + u) 




Ei{9) 
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By the corollary |4.2| and the remark 2.8| we get: 

T2{e) Id 



12 '^S 





Q{6) + Me{e), 



WhercS 



Q{6) :-- 



Xl{0)J{K{e)) 

{j-\K{e))Y{e)yj{K{e)) o 











/« 



and M(.{9) is obtained from Me{9), defined at (4.3), by removing the last n 
rows. So, we have 



(1) 



(4.9) 



T2{e + io) 



Ei{e) 







v-^{e + oj)Q{e + uj) Jj 



Ei{e) + 



+ 



Fj 
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Ei{e) + Me{e)Ei{e), 



(2) 



(3) 



where we used notations: 



v-\e) :- 



Id 



Q{0) 



xi{e)j{K{e)) 
{j-\K{e))Y{e)yj{K{e)) 

Note that, by (4.4) the term (3) in the right hand side of (4.9) is bounded 
by the error i.e.. 



where cq depends on C4 from (4.4) and ||-E'i(^)||J_l Considering (4.5) and an 

elementary computation shows that 

(4.10) 

Fi{d)j-^K{e))Y{e) - Xv{d + uj)Si{e) - j-^{K{e + u;))Y{e + u)-Zv{e + u)A{e) 
F2{e)j-\K{e))Y{e) - XN{e + io)Si{e) - ZN{e + cj)a{9) 



Ei{6) 



substituting (4.10) in the term (1) of left hand side of ( |4.9[ ), we get that 
term (1) is equal to 

xir{e + uj)j{e + uj)E'^'^ 



(4.11) 



v^ 



+ 0J) 



{j'\K{9 + uj))Y{e + oj)yj{e + w)^}'^ 



where E^^ is the upper block of Ei at ( |4.10 ). The definition of Si{6), 
see (2.30), and assumption (2.15) easily show that the lower block in the 



Q{6) is defined at (|2.25[), just to make it easier to follow the calculations we restate it 



again 
5 



As we will see ||i5i(S)||p contains terms that are not bounded by the error, so we do 
not get quadratic bound by the error as in the symplectic case and the constant depends 
on ||_Ei(6')||p also. 



KAM FOR PRESYMPLECTIC MAPS 



23 



equation (4.11) is identically zero. The upper block of the equation (4.11) 



is equal to the following term 

(4.12) (P{e) - i;{9) -Xl{9 + uj)J{e + uj)Xv{e + uj)+ 

- xi{e + uj)j{e + uj)Zv{e + uj)A{e), 

where 

ct>[e) = {F,i9)Xvi9)yipi9)Fiie)J-\K{9))Y{e), 
with ip{9) = J{K{e + w)) - J{f{K{9)) and 

^P{9) = [Fi{9)Xv{9) - Xv{e)VJ{9 + u){Fi{9)J-\K{e))Y{9)). 

Both ^{9) and Fi{9)Xv {9) — Xy (9) are controlled by the error term. This 



fact and Lemma 2.6 show that (4.12) is controlled by ||eo(^)||p. Finally the 

r Q (n\ "I 

term (2) in the left hand side of (4.9) is equal to ^ by definition. 



Since this term is not controlled by the error, we su 







atract it from M ^{9 + 



oj)E{9) to define B{9), then we get the bound (|48j). We move 52(61) to the 
coefficients matrix add it to Si{6). 

U 



Remark 4.4. The details to reach expression (4.12) are as follows: 



Xl{9 + uj)J{9 + u:)J{9 + uj)[Fi{9)J-^{K{9))Y{9) - J-\K{9 + u;))Y{e + w) 
= - [F^9)Xvi9) - Xv{9 + uj)yj{9 + uj)Fi{9)J-\K{e))Y{e) + 



{F^{e)Xv{9)y {J{K{9 + u) - J{f{K{9)) Fi{9)J~\K{9))Y{e) + 



{F^{e)Xv{e)yj{f{K{9))F^{9)j-'{K{e))Y{e) + 



(1) 



Xv{9 + u})J{9 + u)J^^{K{9 + uj))Y{9 + w) . 



(2) 



But we have: 



(1) = n{Fi{9)J-\K{9))Y{9),Fi{9)Xv{9)) = n{J-\K{9))Y{9),Xv{e)) 

(2) = n{J-\K{9 + uj))Y{9 + uj),Xv{9 + w)) = -/, 



so (4.12) follows. 
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We will see that the terms B{e)C{0), Me{e)eo{e) and Me{e){^ 



have a quadratic dependence on the error ||eo 



P' 



^|A=Ao)'^0 

and hence can be con- 



trolled. If we omit these terms from (4.6) we obtain the linear system: 
(4.13) 



Id S{9) 
Id 
A{e) In 



^{e)-^{d + uj) = Ro{e), 



where 

MO) = -V'\9)Q{e)eo{e) - A{e)eo. 
This linear system can be solved using Proposition |2.4[ as we show next: 



Proposition 4.5. Assume that all hypothesis of Lemma 4-3 hold. Then 
there exists a wrapping ^{9), analytic on Up-^25 o-nd a vector eq S M'^" such 



that (4.13) holds for ^{9) andsQ. Moreover, there exits cs and cg depending 

such that 



\p, \\N\\p , 






onn, d, p, r, |/aoIc2,b> II^^o 

(4.14) m\p-25 < csT'S-'^ 



A = An 



|eo| 



(4.15) 



kol < C9|avg(Ao) ^|||eo| 



Proof. Since the proof goes through as in the symplectic case, to avoid un- 
necessary details, we give a short sketch of the proof and refer to [7] for more 
details. Let 



MO) 




m 




so (4.13) becomes 



(4.16) 



f UO) -U0 + ^) = MO) - S{9)^y{9) 

^y{9)-^yi9 + u;) = Ryi9) 
[ UO) - UO + oj) = RM - A{9)iy{9) 



Using the non-degeneracy of the pair {fx,Kx) at A = 0, we can determine 
(sq^ , ...jEq'^) in such way that avg(i?y) = 0. Then we can apply Proposi- 
tion 2.4 to solve the second equation in (4.16) finding a unique zero average 
solution £,y{9). After determining (,y{9) one can choose the remaining com- 
ponents of eo so that 

a-vg{Rx - SCy) 



avg(i?, - A^y) = 0. 
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Applying again Proposition 2.4, we solve the first and last equation of (4.16) 



obtaining unique zero average solutions S,x{d) and ^z{0)- Proposition 2.4 
shows that these solutions satisfy the following estimates: 



UyWp-S < c'^-'6-''\\R 



ywp 



Ux\\p-25 < c""f ^6 "WRx - S^yWp^s 
Uz\\p-25 < c"'j-'S-''\\R, - A^y\\p_S 



The proof of the estimates (4.15) and (4.14) follow just like in the symplectic 
case (see [?)). D 



Corollary 4.6. Assume all the hypotheses of the proposition (4.5) hold, 
then 

(4.17) l|Ao||p-25 < c7-'5-2-||eo||p 



(4.18) \\DGiKo,Xo)\^Ao{e),eo) + eo\\p-2S < 0127"'^^'"+'^ I 

where Ao{0) = M'^ {6)^(6). 



eollp, 



Proof. The estimates (4.17) are immediate consequences of the proposition 



(4.5) and the Cauchy integral formula. Replacing the solution given by 



Proposition |4.5| into the linearized equation (3.8) we find that: 

DG{Ko,Xo)\{Ao{e),eo) + eo{e) = 

' dfx 



M{e + Lo) B{e)i{9) + Mei9)eo{9) + Me{9) 



d\ 



£0 



A=Ao 



Now ( |4.18D follows from (|£7j) (|£8j) (|414j) and (|4T5j). This estabhshes that 
indeed, we have obtained an approximate solution of the linearized equation 
( |3^ . D 

Remark 4.7. One of the concerns in the KAM results of the type we are 
presenting here is how many modifying parameter are needed. A very lucid 
discussion regarding this matter can be found in [25]. A discussion of the 
dimension of the space of parameters in the degenerate cases, can be found 
in [11|. We note that comparing [l\ Proposition 8] and Proposition 4.5, one 



sees, that if the family fx consists of exact presymplectic diffeomorphisms, 
then the dimension of parameter space can be reduced by d. Furthermore, if 
the initial torus satisfies the KolmogorovM non-degeneracy condition il4j i.e. 

It is also known as twist condition 
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if SiYg{S{9)) is non-singular where S{6) is defined in the Lemma 4.3, then 
the dimension of parameter space can be reduced by d again. The reason is 
that we can choose the averages of the tori as parameters. 

In particular, having both families of exact presymplectic mappings and 
Kolmogorov non-degeneracy condition, it will be enough to consider the 
parameter space to be n dimensional, see the Vanishing Lemma also. 

5. Estimates for the improved step 



In the previous section, we have shown that the linearized equation (3.8) 
admits approximate solution in all smaller analyticity domains. The esti- 
mates blow up if the analyticity loss vanishes, the good point is that they 
blow up not worse than a power. 

The goal of this section is to show that if || Ao||p-5 is sufficiently small, the 
new torus Ki{9) = Kq{6) + Aq(6) has an error in the invariance equation 
which is quadratically small with respect to the original one (in the smaller 
domain) . 

Lemma 5.1. Assume 

[Kq + Ao)(C/p_5, Ao + So) C Domain{f), 



where f is defined in (2.3), then 

(5.1) II/ao+.o ° (^"o + Ao) - (Ko + Ao) o T^\l_s < n-^5-^^\\eo\\p, 

where c now involves WfWc^B ^^ ^^^^ ^-^ previous quantities. Furthermore, 
the pair (fx, Ki) is non- degenerate at X = Xq + eo, in the sense of definition 



Note that the linear equation admits estimates for A in any domain Up-s 
for any 5 > 0. If the 5 is very small, the estimates blow up. So that if the 
loss of domain 5 is too small compared with ||eo||p. So that the estimates 
on the step require some restrictions on the of the loss of domain 5 allowed. 



Given the estimates on A,eo obtained in Corollary 4.2, we see that the 
requirement on the composition is implied by 

(5.2) c7-'5-(''^+^)||eo||p<7? 

where r] is smaller than the distance of K[Up) to the complement of the 
domain of /. 
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Proof. This is just a simple consequence of the obvious identity obtained by 
adding and subtracting some terms: 



fXo+eo {Ko + Ao) - (i^o + Ao) o T^ 






iKo)eo-Dh,{Ko)Ao 



A=Ao 



(1) 



+ fx,{Ko) -KooT^ + Dh,{Ko)Ao - Aq o T^ 



dfx 
dX 



{Ko)eo 



A=Ao 



(2) 

The term (1) can be estimated by Taylor theorem, so we have: 

dfx 



\\fx,+eoiKo + Ao)-fx,iKo) 



dX 



< 



1 



A=Ao 
1 



{Ko)eo-Dfx,{Ko)Ao\\p-& 



c^B(l|Ao||p_5 + |eor)<C; 



„-2 r-4o-| 

c2,b7 I 



eol 



The term (2) is exactly the left hand side of (4.18), so by rearranging the 



constant we get estimate (4.18). Non-degeneracy of the pair {fx,Ki) comes 



from the estimates (4.17), (4.15) and the fact that non-degeneracy is an 
open condition. D 

6. Iteration of the Newton method and convergence 
We shall now perform our modified Newton method, starting with /aq. 



Kq, lo and po satisfying the hypotheses of Theorem 3.6, and applying at 



each step the results of Section |4J We will see that if we choose ||eo||po small 
enough we will be able to proceed with the iteration so that the equation 

(6.1) fxiKi9)) = Ki9 + u;) 

has a convergent sequence of approximate solutions 

{Ko,Xo),{Ki,Xi),{K2,X2),... 

defined on domains 

C/po D f/pi D C/p2 D • • • 

with limit an exact solution {Koo,Xoo), defined on a domain Up^. 

Starting with the approximate solution {Ko,Xo), assume that we have 
already found the term {K^, Xm) in this sequence. The next term will take 
the form: 



Km — Km-i + Am-l(^), A^ — A^-l + Em-l {m > 1), 
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with {Am-i{0),em-i) an approximate solution of the Unear equation 

(6-2) DG{Km-l, Xm-l)\{A,r,-i(e),e^-i) = -Cm-l, 

where e^-i := G{Km-i,Xm-i) 



The fohowing lemmas are simply restating the lemma 5.1 for a general 
step. 

Lemma 6.1. Assume that {Km-i, ^m-i) is a non-degenerate (Definition 
|5.5|) approximate solution of (|6.1|) such that 



(6.3) 



rm-i ■= \\Km-i - -K"o||p™_i < r. 



If ||em-i||p„_i is small enough so that Proposition 4-5 applies, then for any 



< Sm-i < Pm-i/3 there exist a function Am-i{6) G 'Ppm-i-ss^-i ^"■^ 
Em-i e M^oi+n^ such that 



Am-1 



'^m-l[0)\\pm-i-25m-i < Cm-ll (^m-^i ||em_l ||p,„_i 

(6.4) \\DA^_i{t 

km-il < Cm-l|(aVg(Am_i)"^|||em_i||p^_j 



-2r-2(f7+l)|| II 



where Cm~i is a constant depending on n, d, p, r, \f\^_^\c'^ s^, \\DKm~i\\p, 



ll^^'fc-illp and 
Moreover if 

(6.5) 



dX 



A — Ati 



1^m-l < Cm-l7 "m-l \\^m-l\\pm-i 



setting Km = Km-i+Am-i, \m = Am-i+em-1- then, em{6) = G{KmAm){0) 
the error function of the improved solutions satisfies 



(6.6) 



-mllp™, ^ f-m-l 7 Wm ll'^m-l llp„_i 



Lemma 6.2. Under the same assumptions as in Lemma 6.1 



one can im- 



prove the constant Cm-i such that (|6.5|) holds and if 

(6.7) 

then 



-2s:-(fT+l)|| II ^ 1 

Cm-l7 "m-l \\^m-l\\pm-l — r. 



(i) // {irDiKm-iy T^DiKm~i is invertihle with inverse N^^i, then the 
matrix {TiDiKmY -KDiKm is also invertihle with inverse N^ satisfy- 
ing 



(6. 



liV, 



m\\pm 



<r WAT II _L ~2r-(o"+l)|| II 

^ IMVm-l||p„_i + Cm-l7 "m-1 ll^m-l ||p™,_i ; 



KAM FOR PRESYMPLECTIC MAPS 



29 



(ii) // Vm-i is invertible then Vm is invertible and the inverse satisfies 

equation ( |6.8[ ) with N replaced by V^^ ; 
(iii) //avg(Am-i) is invertible then avg(Am) is invertible and the inverse 



satisfies equation (6.8) with N replaced by avg(A) 



(1) (iv) The assumption (5.2) ensuring that the range of {Km-i+^m~i, ^m-i + 
£m-i) is inside of the domain of f. 



The most important point is that the constants Cm depend only on n, d, 



\DKm-i\\p, \\Nk-i\\p and 



dX 



A— A„ 



Hence, when 



we show that the K does not leave a neighborhood, then, the constants are 
uniform. 

The convergence of the modified Newton method described above is very 
standard in KAM theory. Indeed, it has sometimes been formulated as an 
implicit function theorem. Among the many versions of implicit function 
theorems, the one of I31j is the closest to the problem here. For the sake of 
completeness, we indicate the main points of the iteration following closely 
[7 16 and refer to those papers for more details. One of the main issues 
to watch out is that the non-degeneracy conditions do not deteriorate much 



along the iteration and that the assumption (5.2), which ensures that we 
can define the composition, remains valid. 

We start by making the choice of the analyticity loss: 



Pm-l 



2-('"~i)(5n. 



The most subtle point is to show that the conditions (6.3) and (6.7) are 



always satisfied. The first one is to guarantee that the new torus always stays 
in the domain of the / and the second one is to insure the non-degeneracy 
condition during the iteration. 

The constant Cm depends on the quantities a, n, d, r, which do not change 
during the iteration. It also depends on the pm < po and the following 
quantities 



I/a„ 



|C2,B^ 



\j^J^m\\pjn, W^mWpmi 



' dX 



I A — A77 



(K. 



mJWPm 



, I avg(A^ 



This dependence is polynomial. By similar calculation as follows, can be 
shown that there exist constant c such that Cm ^ c for m, > 0, see [tI Lemma 
13]. The main point is that we do not get far away from initial torus. Denote 
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Em = 1 1 Cm lip™, with the choice of the domain loss, we obtain : 

(6.9) 

6™ < c7-^(2-(™-i),5o)-^'^e^_i < (c^-4)(i+2)(2-K-i)-2(™-2)^(i+2))-4. 

< (c7-^<5o-"")''"-^24-(2'"-^m)^2'" < (c7-4^-4a24^,p)2--l2-4<x(m-l)^^^ 

where we have used that 



-4o- 4 

^m-2 



m— 1 



2°(m - 1) + 2(m - 2) + • • • + 2"^-2 = 2"^-^ ^ s2-" < 2^ 



m. 



s=l 



One sees that if ||eo||po satisfies the assumption (3.3), the condition (6.7) is 



always satisfied. It remains to show that (6.3) is also satisfied. We denote 



K = cj ^(5q '^2^'^eo. Now, the first estimate in (6.4), estimate (6.9) and the 

definition of r^ gives us: 

(6.10) 



m— 1 



fm _; ^m—1 ~T~ Cm'J 0,^— 1 ll^"^~l llpn 



^ -2 -2o- , -2 -2a 



,<■■■< c7-Vo-2-eo + C7-' J] <^-'"e,- 

m— 1 



c7^'^o"'"eo 



[l + -2- £ 2- A = c7-Vo--eo (l + .^) . 



Again having the assumption (3.3) and calculations (6.10) show that the 



(6.3) is always satisfied. 



7. Local Uniqueness 



The proof of Theorem 3.7 follows exactly the same pattern as the proof 
of uniqueness for the symplectic case given in |7] , so we will not reproduce 
here the same computations. We limit ourselves to some comments and a 
sketch of the proof, which takes advantage of the fact that in Proposition 



2.4 two different solutions of (2.8) differ by their average. In our situation. 



one can transfer this difference of averages of two solution to a difference of 
the phase between them. 

Let {Ki, Ai) and {K2, A2) be two solutions as in the statement of Theorem 



3.7 From Taylor's Theorem we have: 

(7.1) Dfx,iK{6)){K2-K,)+R{Ki,K2) = 
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where 

(7.2) 



\\R{K,,K2)\\p<c\\K2-Ki 



Applying the change of variable MS, = {Ki — K2), where M is given by 
(|2.17) and replacing K by Ki, the linear equation (7.1) is transformed to 



(7.3) 
(7.4) 
(7.5) 
where 



UO) - UO + ^) = {R{Ki,K2)h - Si{6)iy{9) 
Cy{0)-CyiO + i^) = iR{Ki,K2))y 

UO) -uo + oj) = iRiKi,K2)), - A{e)Sy{e) 

RiKi,K2) = -M-\e + u)R{Ki,K2). 



Using Proposition 2.4 and (7.2), it follows from (7.4) that: 
(7.6) Uy\\p-25 < n-'5-^''\\R{KuK2\\p 

where 



tid) = CyiO)-e^yg{Q 



On the other hand, the average of the right hand sides of (7.3) and (7.5) are 



zero, so the assumption that (see Theorem 3.7) has rank d together with 



the estimates (7.2) and (7.6) give us 



\m,Cy,Cz) - (avg(e.),0,avg(e,))||^_2,5 < cj-^S-^^\\K2 - Ki\\l 
Similarly to |7|, this leads to the following lemma. 

Lemma 7.1. There exists a constant c depending ond, n, p, |J|b,., ||-f^i||c2,p 
such that if 

c\\K2-Ki\\p<l, 

then one can find r S R +" with \t\ < \\K2 — i^ilL such that 



avg 






[K2 



^illp 
i^i]^ = 0, 



where Ti and T3 are defined by (2.29) after replacing K by Ki. Therefore, 
for any < 5 < p/2, we have 

\\Ki o Tr, - K2\\p-25 < CJ-^d-^^'WKi - K2\\p, 

for a constant c depending on the same parameters as c and also on G^^. 

We then replace K2 by K2oTt-^ and repeat the iteration, which is possible 
since K2 o T^-^ is also an invariant torus and the constants c and c do not 
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depend on K2. In this way we produce a convergent sequence of phases 
n, r2, . . . , Tm, • • • such that the hmit Too satisfies: 

\\K2 0Tr^-K4p/2 = Q, 



therefore completing the proof of Theorem 3.7 



A different proof which does not require iteration (but requires setting a 
normahzation condition) appears in (sl. 



8. Comparison between Poisson and Presymplectic cases 

As we mentioned in the introduction, our results are not applicable to 
the Poisson dynamical systems. In this section, we only present a sim- 
ple example to justify our statement about the difference between Poisson 
and presymplectic diffeomorphisms. Consider the presymplectic structure 
il = dx A dy on M := T*T x T with standard coordinates {x,y,z). The 



corresponding compatible Poisson structure is 11 



d 



A 



d 



dx dy 



For the diffeo- 



morphism / : M — )• M to persevere the Poisson bivector 11, it has to satisfy 
following condition: 



^.1) 



(Df) 






-1 







-1 


1 





{Dfy = 


1 




















which means / will have the form 

f{x,y,z) = {fi{x,y,z),f2{x,y,z),h{z)), 



and ^^ 
^^°- dx dy 



ddx ~ ^^ ^^ note that \S.1\ clearly shows the symplec- 



tic leafs of M are invariant under /, which is a simple fact from Poisson 
geometry. For / to be presymplectic diffeomorphisms of M, it has to satisfy 



{Dfy 

this leads / to have the form 




1 




-1 







{Df) 






-1 





1 


















f{x,y,z) 
dfi a/2 a/i a/2 



{fiix,y),f2ix,y)j3{x,y,z)), 



and again ^q^^q^ — ^ ^^ = 1. In general, there is no canonical way to get 
a symplectic foliation for a presymplectic manifold. Even if we consider the 



symplectic foliation arising from Poisson structure, one can see from (8.2) 



that presymplectic diffeomorphisms, in general, do not preserve symplectic 
leafs. 
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